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PACS 78 . 47 . -p - Spectroscopy of solid state dynamics 

PACS 74.50.+r - Tunneling phenomena; point contacts, weak links, Josephson effects 
PACS 42 . 50 . Pq - Cavity quantum electrodynamics; micromasers 

Abstract. - We show that the nonlinear response of a driven circuit quantum electrodynamics 
setup displays antiresonant multiphoton transitions, as recently observed in a transmon qubit 
device. By including photon leaking, we explain the lineshape by a perturbative and a semiclas- 
sical analysis. We derive a bistable semiclassical quasienergy surface whose lowest quasienergy 
eigenstate is squeezed, allowing to define a squeezing-dependent local effective temperature. We 
study the escape dynamics out of the metastable state and find signatures of dynamical tunneling, 
similar as for the quantum Duffing oscillator. 
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Introduction. — One of the nontrivial fundamental 
models of quantum physics is the Jaynes-Cummings (JC) 
model pp. It was introduced to describe the interaction 
of a two-level atom and a single quantized electromag- 
netic field mode. Being sufficiently simple, its dynam- 
ics is very rich though, including Rabi oscillations, col- 
lapse and revival phenomena, squeezing, entanglement, 
Schrodinger cat and Fock states, and photon antibunch- 
ing [2]. Beyond quantum optical set-ups, it is applicable 
to many situations of nanocircuit quantum electrodynam- 
ics (QED), such as Cooper pair boxes [3], superconducting 
flux qubits [4] , Josephson junctions [5] , and semiconduc- 
tor quantum dots In particular, the latter setups allow 
to explore the regime of strong coupling and nonlinear re- 
sponse. 

Recently, unique nonlinear features have been detected 
in the transmitted heterodyne signal of a superconduct- 
ing transmon qubit device [7]. For weak driving, the two 
well-known vacuum Rabi resonances reflect transitions be- 
tween the JC groundstate and the first/the second ex- 
cited state of the undriven JC spectrum. Their difference 
in energy is 2hg, where g is the interaction strength of 
the qubit and the harmonic mode. For increasing driv- 
ing, each vacuum Rabi peak supcrsplits into additional 
(anti-)resonances with the characteristic ^/n spacing. The 
measurements have been corroborated with accurate nu- 
merical simulations based on a Markovian master equa- 
tion of Lindblad form [7]. For the particular case of 



the one-photon resonance of the JC groundstate and the 
antisymmetric superposition of transmon and photon ex- 
citation, an effective two-level model composed of both 
transmon and cavity degrees of freedom has been estab- 
lished. The ensuing Bloch equations have been solved an- 
alytically, yielding a closed expression for the heterodyne 
amplitude. For higher multiphoton resonances, the nu- 
merical solution of the master equation was compared to 
the experimental data, with superb agreement. However, 
the underlying nontrivial competition of driving and dis- 
sipation has not been discussed. 

In this Letter, we provide a complete physical picture for 
the dissipativc dynamics of the Purcell-limitcd nonlinear 
response of the driven JC model in terms of quantum mul- 
tiphoton (anti-)resonances. The underlying competition of 
driving and dissipation is revealed by perturbative argu- 
ments in the rotating frame in presence of photon leaking. 
We show explicitly that the lineshape of the multipho- 
ton resonance in general is determined by the ratio of the 
Rabi frequency and the dissipation strength, allowing for 
a direct experimental control. In addition to the general 
analysis of the lincshapes of the multiphoton resonances, 
which were discussed in Ref. [7] based on the numerical so- 
lution of the master equation only, we consider the semi- 
classical limit of the driven JC model. In the regime of 
strong driving, which cannot be accessed perturbatively, 
we derive a novel semiclassical quasienergy surface which 
is bistable. We show that its lowest quasienergy eigen- 
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state is an amplitude squeezed state and displays large 
out-of-phase oscillations. We consider the dissipative dy- 
namics on the quasienergy surface at zero temperature. 
It is known that dissipation in driven nonlinear oscilla- 
tors introduces nontrivial effects, such as quantum activa- 
tion [8"HTU] and dynamic resonant tunneling [TTMl3j . We 
show that for the case of the dissipative driven JC model, 
the lowest quasienergy eigenstatc is significantly popu- 
lated at a multiphoton (anti)resonance, and is mctastable 
away from resonance. The scmiclassical analysis allows for 
an intuitive picture of dynamic bistability and dynamic 
resonant tunneling for the driven dissipative JC model, 
which has not been developed before. Furthermore, we 
reveal deep analogies with the quantum Duffing oscilla- 
tor 0HUHI6]. 

The model. — We start from a harmonic oscilla- 
tor with frequency w r which is resonantly coupled with 
strength g to a qubit of the same frequency ui T and which 
is driven with frequency w GX and field strength /. In the 
frame rotating with w cx and for 5lo = ui x — lo cx , g, f <C w r , 
we perform a rotating-wave approximation and obtain the 
Hamiltonian of the driven JC model (h = ks = 1) 

H = Su Ua + y + 5) +f (« tfT - + a °+) + { ( flt + a ) > 

(1) 

with a± = <j x ± itjy. Here, aj arc the Pauli matri- 
ces. The undriven JC model has the quasienergies Sq = 
—5uj/2 , £ n ,± = (n — 1/2)Suj ± g\/n and the quasienergy 
states \4> ) = |0, 5 ), \<j> n ±) = {\n-l,e)±\n,g))/V2. We 
will refer to latter as n-photon (dressed) states with two 
spin directions ±. For / 7^ 0, avoided crossings of the 
quasienergy levels arise, which correspond to ./V-photon 
transitions at 5uj = ±g/\/N. To have well separated 
resonances, we consider the regime g /. Around the 
resonance, the zero-photon state |^>o) and the iV-photon 
dressed state \4>n±) display Rabi oscillations with the Rabi 
frequency ilpj oc f N /g 1 ^^ 1 . 

Due to the suppression of 1//— noise, transmon qubits 
are Purccll-limitcd in the resonant regime considered 
here [17] . For low temperatures T <C w C x; we incorpo- 
rate Purcell dissipation by means of the simple Lindblad 
master equation 

p = -i[H, p} + ^ ([ap, a f ] + [a, pa)]) , (2) 

for the reduced density operator p of the coupled qubit- 
oscillator system. Then, only photon leaking from the 
system into the bath is possible. Hence, in absence of mul- 
tiphoton transitions, |0o) is dominantly populated in the 
stationary state. In contrast, at a multiphoton resonance, 
the stationary state is generated by coherent driving to the 
-/V-photon state and a subsequent relaxation via all the in- 
termediate n-photon states (n < N) to the 0-photon state 
due to photon leaking. Eventually, this nontrivial inter- 
play generates a stationary mixture of all rt-photon states. 




Fig. 1: Amplitude A as a function of the driving frequency ui ex 
for / = 0.0006cj r (a), / = 0.0015oj r (b) and / = 0.0025cj r (c). 
Insets: zooms to the corresponding (anti-)resonances marked 
by the dashed lines. Moreover, g = 0.026oj r ,K = 5 x 10 u> r . 
A realistic value for o; r /(27r) is 7 GHz [71IT7]. 



In particular, each different choice of model parameters in- 
fluences the intermediate dissipative transitions, which, in 
turn, determine the asymptotic populations of the station- 
ary state. 

We are interested in the nonlinear response character- 
ized by the steady-state intracavity field (a) = tr(ga) = 
0ap a f3a, where a,, ft refer to the quasienergy 
eigenstates basis. We discuss the case 5ui > 0, the opposite 
follows from \<f> n ±) — > \4> n -^), f ^ — f and (p — > —<p + n. 
The modulus A is related to the transmitted amplitude 
Atr oc A and intensity 7t r oc A 2 of the microwave input 
signal at frequency w ox [7]. In the rotating frame, (a) < 
(ip = 7r) corresponds to an oscillation out of phase with 
respect to the drive. 

Numerical results. — The nonlinear response, in the 
first instance obtained from a numerical solution of Eq. 
([2]), is shown in Fig. [1] for the parameters in the regime 
of Ref. [7]. The drive induces a splitting of the vacuum 
Rabi resonance and produces two families of peaks which 
are symmetric with respect to iu ex = w r and which are 
associated to the ± quasienergy states. In each family, 
(anti)resonances occur which correspond to multiphoton 
transitions and which are associated to the avoided cross- 
ings of quasienergy levels, see Fig. [5b. We note that no 
antiresonances occur in the photon number (a) a) [18] . 

For weak driving (Fig. QJi), only the 1-photon antires- 
onance is well pronounced. It can be described [7J[T3] by 
a model involving the 0- and the 1-photon state. The 
generic mechanism is that, at resonance, they are equally 
populated and oscillate with opposite phase yielding zero 
response. Slightly away from resonance, one of the two 
states is more populated and a finite response arises. Far 
away from the resonance, the response again approaches 
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Fig. 2: Nonlinear response of the driven JC model: (a) ampli- 
tude, (b) phase, (c) quasienergies, and (d) population of the 
lowest quasienergy state for g = 0.026cj r , / = 0.004a; r , re — 
10 -4 cj r - Dashed-dotted red line in (a,b): lowest-order result 
for non-resonant approximation. Dashed orange line in (a): 
A* = |(^* |a|^*)|; and in (c): semiclassical result eq. ([6]) for 
the lowest quasienergy E* . 



zero. Overall, the lineshape of an antiresonance arises. 
The antiresonance around u) ex ps 0.98w r corresponds to a 
2-photon process (Fig. ^,). This feature also follows from 
a two state description. The shape is different from the 1- 
photon antiresonance, due to the background contribution 
of nonrcsonant 1-photon mixing processes |16j . 

For increasing driving, unexpected features arise. The 
multiphoton antiresonances turn into resonances, see 
Fig. Hp and c. Associated to the behavior of A are jumps 
in the phase ip, see Fig. [5b. This already suggests that 
two quasienergy states - one oscillating in and one out of 
phase - arc alternately populated. Interestingly, the pop- 
ulation (Fig. [2Ji) of the state \tp*) with lowest quasienergy 
shows peaks at the resonance frequencies. In fact, as we 
will show below, is localized in the bottom of a well 
of a bistable quasienergy surface. \ip*} is metastablc since 
the bath induces transitions to higher quasienergies states 
and an escape is always possible, even at zero tempera- 
ture. This feature has also been reported for the quantum 
Duffing oscillator [TTHI3I- 

Perturbative approach. — These observations are 
further substantiated by perturbative arguments, which 
show how the asymptotic populations of the stationary 
state are determined by the various rates of the interme- 
diate dissipative transition processes. Here, we assume 
that the dressed states \4> n ±) (n ^ N) are approximate 
eigenstates. In the next section we will show that this is 
true only for very weak driving fields. 

Out of resonance, pqq — 1 yielding (a) ~ —f[(5u) + 
g) _1 + (5uj — g)^ 1 ]/4 up to first order in /. As follows 
from figs. [H a) and b) (dashed-dotted red lines), the lowest 
order response A$ coincides with the exact one away from 
resonance. Moreover, ipo = n for Slo < g and ipo = for 
5uj > g. 



At the iV-photon-resonance, the system can leave \cj>o} 
by tunneling to \4>n—) and the bath induces decays from 
\4> n -} to |(/)„_i_) along the ladder N — > N — 1 — > ... — >• 
1 — > 0. The decay rates (in secular approximation) are 
readily evaluated from Eq. as C n _/ n _i\_ = (y 7 ?! + 
yjn— l) 2 re/4 for i! ^ 1 and £i_,o = k/2. Hence, the 
rate from the 1-photon to the 0-photon state is smallest. 
Note that the probability of a decay to a state with oppo- 
site dressed spin is small, i.e., £( n -i)-. «+/£(«-!)+■«+ — 
l/[16n (n — 1/2)]. Hence, the population of the 1-photon 
state is always larger than those of the n-photon states. In 
addition, depending on the ratio k/Qn, qualitatively dif- 
ferent stationary populations arise from a competition be- 
tween tunneling from | </>o) to the top of the ladder, |0jv_), 
and relaxation down the ladder. For re 3> 51 at, the pop- 
ulation of \(f>o) is ~ 1, because damping is more efficient 
than tunneling. Dissipation then completely washes out 
the resonance, and the response is the same as that off res- 
onance and thus is in phase with the drive. For re ~ f2jv, 
a small population p\\ emerges, contributing p\\<X\\ with 

= 2( El _-e^) 3+ 4^ + 4( £l l f - £0 ) < sinCG ' in the P Cr " 

turbative regime and for N < 6, S\- < £2-,£o- This 
leads to a reduced response, forming an antiresonance, 
see, e.g., the 2-photon antiresonance shown in Fig. [T]i. 
For re <SC fijv, tunneling is faster than relaxation and the 
population of |0i_) becomes the largest. Then, the contri- 
bution pnan < dominates, leading to an out-of-phase 
oscillation. To be specific, for the 2-photon resonance of 
Fig.Q] k/Q 2 = ng(V2f 2 ) = 1.8 (a), 0.29 (b) and 0.1 (c). 

Semiclassical approach. — For increasing (but still 
weak) driving , the response is qualitatively similar, al- 
though the perturbative approach becomes inadequate. In 
fact, as shown in Fig.jSJi, at resonance, there is a large pop- 
ulation P* of the lowest quasienergy state \ip*). As follows 
from Fig.^, \ip*) ^ \<fii-)- In order to account for the im- 
portance of \ip*}, we perform next a semiclassical analysis. 
We switch to the dressed atom picture, by means of the 
unitary transformation 

R = cxp [ =[qt<7_ — a<r + ] | , (3) 

^8^0 + ^/2 + 1/2 + 7 

mapping the JC eigenstates into product states \na) (with 
a = g,e). Under this transformation, the undriven JC 
Hamiltonian becomes diagonal in the spin degrees of free- 
dom and reads 



H = \8u\la*a + -£ + -j+g<r z ^a!a+-£ + - (4) 

while a = i? f ai? = a + O^ 1 / 2 ). This expression follows 
by expanding the matrix elements of the creation opera- 
tor in the dressed Fock basis with respect to n. Higher 
order terms include spin flipping operators. Next, we in- 
troduce the canonical variables X = ^/A/2(o^ + a) and 
V = iy/\/2(a ! — a), where A = Suj 2 /g 2 is a dimensionless 
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Fig. 3: Contour plot of the quasienergy surfaces Q(X,V) for 
the dressed spin orientation g (a z — —1) (a) and e(a z — +1) 
(b) for / = 0.004oj r and g = 0.026aj r , yielding f/g = 0.154 and 
Teft = 0.65. The solid black lines in (a) indicate quasiclassical 
orbits and include also the separatrix between the two domains 
of attraction. 



parameter playing the role of H. Eventually, neglecting all 
higher order terms , we obtain the transformed Hamilto- 
nian H ~ 5ujX~ 1 Q{X, V) with 



X 2 V 2 



f 

— + — + -i^x 



X 2 
2~ 



■p2 

2 



(5) 



It can be interpreted as a quasienergy surface in phase 
space and is visualised in Fig. [3J For a z = —1 (Fig. [3] a), 
an internal dome appears around X = 0, V = 0, such that 
two quasiclassical orbits encircle the inner maximum, one 
on the external and one on the internal domain. They are 
degenerate with respect to their quasicnergies and thereby 
define a dynamical bistability. Both domains are sepa- 
rated by the separatrix shown in Fig. [3] a as the oriented 
black solid line. The overall shape of the quasienergy sur- 
face resembles a Mexican hat. The drive induces a finite 
tilt of the surface, generating a quasienergy well around 
the quasienergy minimum. The orbits near the inner max- 
imum correspond to small photon numbers and are not ac- 
curately described by the semiclassical approach. On the 
other hand, the orbits on the outer domain correspond 



to states with large photon numbers. Deep in the quan- 
tum regime, discrete multiphoton transitions correspond 
to tunneling transitions between these quasiclassical or- 
bits. The surface for the opposite dressed spin orientation 
a z = +1 is a less interesting monotonous function, shown 
in Fig. [3Jb. 

Let us next consider the dynamics around the minimum 
located at V m i n = and X min = — (f/g + 1) /\/2. A 
harmonic expansion yields the quasienergy E* for the state 
localized at the bottom of the quasienergy surface as E* = 
5u}\~ 1 [Q(X min , 0) + (A/2)w*], or written in terms of / and 
9 as 

I/O \ 2 

9 

4<5w 



E* 



J - + l\ +^ 




9+r 



(G) 



with the effective frequency uj* = \/ f/{g + /)• This result 
is correct up to 0(\ 2 ) and is shown in Fig. [5p as orange 
dashed line. It almost coincides with the exact result, even 
for small photon numbers N = 1. 

Up to leading order in A, the lowest quasienergy eigen- 
state 

\r) ^R^D(X mm )S(r)\0g) (7) 

is obtained by applying to the vacuum: i) the squeez- 
ing operator S(r) = exp [r(a 2 — a f2 )/2\ with squeeze fac- 
tor r = -(lnw*)/2 = ln[l + g/f}/4, ii) the translation 
D(X m - m ) = exp [iVX m - m / X] to the minimum, and iii) Rr 1 
to return to the bare atom picture. With this, one can 
compute the Fock-state representation Q15] and all expec- 
tation values at leading order. It turns out that has 
sub-Poissonian statistics and shows photon antibunching. 
This is not surprising, in fact D(X m i n )S(r)\0 g) is an am- 
plitude squeezed state for r > 0. 

Note that this state \ip*) exists for any finite driving 
strength (provided that A be much smaller than the bar- 
rier height f/g), but it is not present in the undriven case. 
For this reason it is not possible to describe it in terms 
of simple perturbation theory in the driving. In fact, in 
the limit of vanishing driving, the tilt of the quasienergy 
surface vanishes and no quasienergy well is developed. For 
weak but finite driving, the well is still very shallow in the 
momentum direction, allowing for large momentum fluctu- 
ations of a state confined in it. For increasing driving, the 
well becomes deeper and more symmetric. For this rea- 
son, the squeezing factor decreases for increasing driving. 
Let us furthermore emphasize that \ip*), i.e., the lowest 
quasienergy eigenstate of the driven JC Hamiltonian, can 
not be identified with the squeezed state investigated in 
Ref. [20] (for 5cj = 0). 

The latter is the eigenstate obtained by starting in the 
JC groundstate \4>o) and by switching on adiabatically the 
driving. It has vanishing quasienergy for 5u — > 0. More- 
over, it tends to follow the driving by rotating its spin 
part and squeezing its amplitude fluctuations. Hence, as 
opposed to \ip*}, its squeezing factor increases for increas- 
ing driving. We expect the Wigner representation for this 
state to be confined in an area of size A 2 around the local 
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W ex /(0 r 

Fig. 4: Solid line: Smallest eigenvalue of the Lindblad master 
equation. Dashed red line: the same without dissipative spin 
flips. Dashed-dotted orange line: decay rate for \ip*) — ► \</>q). 
Dotted lines correspond to k/2 and «/4. The parameters are 
g = 0.026cj r , / = 0.004aj r , and k = KT 4 o; r . 



quasienergy maximum. Unfortunately, it cannot be incor- 
porated in our semiclassical description because it corre- 
sponds to a small photon number. In the following dis- 
cussion, we will assume that / <C g, so that perturbation 
theory is still valid close to zero quasienergy (although it is 
not valid close to the quasienergy minimum). In this case, 
the JC groundstate \(f>o) is still an approximate cigcnstate 
of the driven JC Hamiltonian. 

Next, we consider the dissipative semiclassical dynam- 
ics. As will be shown below, a separation of time scales 
exists which defines a fast intrawell and a slow inter- 
well relaxation. Deep in the semiclassical limit, the 
quasienergy states, localized close to the minimum of one 
well, can be obtained as \ip*) = J2 — 1 6 t " R\ip*) / '\fn\ with 
b = a cosh r + sinhr — e r X m i n /-\/2\, In this limit, dis- 
sipative transitions occur only between nearest neighbors, 
with the rates C n -i. n = «ncosh 2 r, C n , n -i = Knsinh 2 r. 
Here, the detailed balance condition is fullfilled. When the 
system is initially in a state with a large photon number, 
it has a large probability to fall in the basin of attrac- 
tion of the quasipotcntial minimum (intrawell relaxation) . 
When also k <C gAw*, detailed balance determines an ef- 
fective Boltzmann distribution P* = p*e~ nt3ot{ , with effec- 
tive inverse temperature /3 e ff = 21ncothr. We emphasize 
that this link between effective temperature and squeezing 
can be generalized to any driven quantum system with a 
smooth quasienergy surface and coupled linearly to a bath 
(e.g., the Duffing oscillator [S UTDT - fTB] or the parametrically 
driven oscillator [3]). It can be easily generalized to fi- 
nite temperatures T > as well. It turns out that the 
zero temperature limit applies when sinh 2 r is much larger 
than the bosonic occupation number n(uj ex /T) of the bath 
at Wex- In the opposite limit, /3 e ff = ^ck/T. Since we in- 
clude here only photon leaking, i.e., uj cx 3> T, the effective 
temperature is still small. 

On the large time scale, the system decays to the 0- 



photon state with a forward rate k+ for the interwell re- 
laxation. From there, it can return to the basin of attrac- 
tion of the minimum by a driving induced transition with 
a backward rate fc_ . The stationary populations of the in- 
trawell states (oscillating out of phase) and the 0— photon 
state (in phase) are determined by the ratio k+/k-. Off 
resonance, photon leaking favors the 0— photon state and 
k + 3> k-. Approaching a resonance, fc_ increases due 
to an increasing contribution due to tunneling to the es- 
cape. When the tunneling contribution becomes compa- 
rable to k + , the response is qualitatively modified yielding 
an (anti-)resonancc. 

To render the qualitative discussion more quantitative, 
we consider the eigenvalues of the Liouville operator which 
governs the Lindblad master equation ([2]). We start deep 
in the semiclassical regime, i.e., for small detuning, where 
uj cx pa ui r . Here, a clear separation of time scales for the 
dissipative dynamics on the bistable quasienergy surface 
occurs. Well-defined energy wells with a large quasienergy 
barrier in between exist for small detuning and allow for 
a clear description in terms of a single relaxation rate (see 
also Ref. [U). I n this regime, we find a single eigenvalue 
r which consists of the sum of fc_ and fc+, is real and 
much smaller than the real parts of all the other eigen- 
values. This is shown in Fig. [4] as black solid line. For 
increasing detuning, i.e., when w cx < 0.987o; r in Fig. [4l 
we enter a regime, where the separation of time scales is 
not so clearly expressed and V becomes comparable to the 
real parts of three more eigenvalues. The latter corre- 
spond to relaxation (one real eigenvalue) and to decoher- 
ence (a complex conjugated pair of eigenvalues) involving 
the pair of states \<j>\+) and \4>o). We do not show them 
in the figure, but instead show the perturbatively deter- 
mined values k/2 (relaxation) and k/4 (decoherence) out 
of resonance as dotted horizontal lines. The peaks in Fig. 
[4] (black solid lines) are due to tunneling from the 0- to N- 
photon state. When the Rabi frequency fi^r is much larger 
than the A— photon dressed state decay rate T^r, the sys- 
tem tunnels back and forth many times between the two 
internal domains of the quasienergy surface, before it de- 
cays to the external domain with rate fc_ ~ Yn /2. In the 
opposite limit Tn ^> Qn , the tunneling rate is reduced due 
to the quasienergy level fluctuations. It becomes equal to 
flfj /r n [22] and coincides approximately with the escape 
rate k-. Away from resonance, r w fc+. 

There are three mechanisms of decay from the 
metastable well: (i) The system can decay directly to the 
0-photon state with the rate shown in Fig. @J dashed- 
dotted orange line, (ii) The system can climb up the 
quasienergy well by quantum activation [jjj. Both asso- 
ciated rates are expected to decrease exponentially, fol- 
lowing tx e~ Ci / Su , with some constants q (the prefactor 
varies smoothly with Soj), which defines the separation of 
time scales, (iii) For very small detuning, the escape oc- 
curs via bath-induced spin flips. In fact, this mechanism 
is suppressed only as a power-law T oc Suj 2 . To separate 
(ii) from (iii), we show T without the bath- induced spin 
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flips in Fig. @] (dashed red line), illustrating that spin flips 
are dominant when the separation of time scales is well 
defined. Since only a few states close to the bottom are 
populated, our solution is stable against a small dephasing 
of the oscillator [S] or an intrinsic spin relaxation that vi- 
olates detailed balance. The induced spin-flip rate would 
be small and remain finite for A — > 0, imposing an upper 
limit to the lifetime of \4>*)- 

Perspectives and conclusions. The transforma- 
tion of the driven JC Hamiltonian to the quasienergy pic- 
ture defined by Eq. ((5]) allows for a quantitative semiclas- 
sical WKB analysis in the spirit of Refs. [8j[9]. We defer 
this to a future publication since it lies beyond the scope 
of this Letter. Our analysis can be easily extended to any 
driven nonlinear oscillator coupled bilinearly to a thermal 
bath. For example, the Duffing oscillator is characterized 
by two classical stable solutions with different amplitudes 
and opposite oscillation phase. The two states can be as- 
sociated with classical orbits in a quasienergy landscape 
which also assumes the form of a tilted Mexican hat [T3] 
like for the JC model. Transitions between the two states 
can be induced by three mechanisms: (i) thermal activa- 
tion "over the barrier" [MlfPo] , where the thermal fluctua- 
tions provide the energy for classical escape, , (ii) quantum 
activation [5j where bath zero-point quantum fluctuations 
induce a diffusive motion on the quasienergy surface over 
the barrier, which even at zero temperature can lead to 
escape, and (iii) quantum tunneling [11H13] . 

For weak driving, the small-oscillation solution can be 
identified with the 0— photon quantum state. Since it is 
favored by photon leaking, it has a low effective temper- 
ature and can be regarded as stable in absence of tunnel- 
ing. However, this stable state is associated to a relative 
quasienergy maximum, in contrast to the case of a static 
bistable potential. At a multiphoton resonance, this state 
becomes mctastablc and generates an (anti-)rcsonance of 
the stationary oscillation and tunneling peaks in the inter- 
well relaxation rate as a function of the detuning. These 
features has been already predicted in Refs. fTTHTB"] . but 
the link to the semiclassical picture was not drawn. 

Note that our results on the resonant tunneling transi- 
tions in interwell dissipativc dynamics is not in contradic- 
tion with the findings of Ref. [5] . There, it has been shown 
that the quantum activation rate is always larger than that 
of dissipative tunneling when they are computed semiclas- 
sically with logaritmic precision. Here and in Refs. [TTHT3] . 
we show that tunneling is the dominating escape mecha- 
nism in the complementary deep quantum regime when it 
is coherent and/or only very few states are confined within 
the quasienergy well (so that logaritmic corrections are rel- 
evant). 

We furthermore note that the dissipative quantum Duff- 
ing oscillator could be realized in a Josephson bifurcation 
amplifier |23j operating in the quantum regime. 

In conclusion, inspired by recent experiments, we have 
explained the nonlinear response of the driven dissipative 



Jaynes-Cummings model. We have predicted the existence 
of a metastable squeezed state in the semiclassical limit 
and drawn a link between effective local temperature and 
the squeezing parameter. We have analyzed the escape 
mechanisms from the metastable states and found reso- 
nant dynamical tunneling. Our analysis reveals generic 
features on the dissipative dynamics of nonlinear driven 
quantum systems for the example of the driven Jaynes- 
Cummings model, realized in terms of a nanocircuit QED 
setup. 
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